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Abstract —ff an interarea oscillatory mode has insufficient 
damping, generator redispatch can he used to improve its 
damping. We explain and apply a new analytic formula for 
the modal sensitivity to rank the best pairs of generators to 
redlspatch. The formula requires some dynamic power system 
data and we show how to obtain that data from synchrophasor 
measurements. The application of the formula to damp interarea 
modes is explained and illustrated with interarea modes of the 
New England 10-generator power system. 

Index Terms —Power system dynamic stability, phasor mea¬ 
surement units, power system control. 

I. Introduction 

Power transmission systems have multiple electromechani¬ 
cal oscillatory modes in which power system areas can swing 
against each other. In large grids, these interarea oscillations 
typically have low frequency in the range 0.1 to 1.0 Hz, 
and can appear for large or unusual power transfers. Poorly 
damped or negatively damped oscillations can become more 
frequent as power systems experience greater variability of 
loading conditions and can lead to equipment damage, mal¬ 
function or blackouts. Practical rules for power system security 
often require sufficient damping of oscillatory modes Cl, 12, 
such as damping ratio of at least 5%, and power transfers on 
tie lines are sometimes limited by oscillations Cl, 0, a, 0. 

There are several approaches to maintaining sufficient 
damping of oscillatory modes, including limiting power trans¬ 
fers a, installing closed loop controls [Tl, and the approach 
of this paper, which is to take operator actions such as redis¬ 
patching generation ©, Q, ®. It is now feasible to monitor 
modal damping and frequency online from synchrophasors 
(also known as PMUs) a, cni- Suppose that a mode with 
insufficient damping ratio is detected. Then what actions 
should be taken to restore the mode damping? 

This paper calculates the best generator pairs to redispatch 
to maintain the mode damping by combining synchrophasor 
and state estimator measurements with a new analytic formula 
for the sensitivity of the mode eigenvalue with respect to 
generator redispatch. This formula, previously thought to be 
unattainable, is derived with a combination of new and old 
methods in im. The length of the derivation (more than 8 
pages) precludes its presentation here. In this paper we state, 
explain, and demonstrate the application of the new formula 
and show how the terms of the formula could be obtained 
from power system measurements. In particular, we propose 
using synchrophasors to measure the terms of the formula that 
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depend on dynamics and using the state estimator to measure 
the terms of the formula that depend on statics. In applying 
a first order sensitivity formula, we assume that the power 
system ambient or transient behavior is dominated by the 
linearized dynamics associated with an asymptotically stable 
operating equilibrium. 

Changes in generator dispatch change the oscillation damp¬ 
ing by exploiting nonlinearity of the power system: changing 
the dispatch changes the operating equilibrium and hence the 
linearization of the power system about that equilibrium that 
determines the oscillatory modes and their damping. This open 
loop approach that applies an operator action after too little 
damping is detected can be contrasted with an approach that 
designs closed loop controls to damp the oscillations preven¬ 
tively. The closed loop control design chooses control gains 
that appear explicitly in the power system Jacobian, whereas 
generator redispatch changes the Jacobian indirectly by chang¬ 
ing the operating point at which the Jacobian is evaluated. 

We now review previous works using generator redispatch to 
damp oscillations; these have considered heuristics, brute force 
computations, and formulas that are difficult to implement 
from measurements. These approaches have established that 
generator redispatch can damp oscillatory modes. 

Fischer and Erlich pioneered heuristics for the redispatch 
in terms of the mode shapes for some simple grid structures 
and for the European grid 0, M- Their heuristics seem 
promising for insights, elaborations and validation, especially 
since there has also been progress in determining the mode 
shape from measurements II3, m, Ga, GSl. 

There are also previous approaches that require a dynamic 
grid model. The effective generator redispatches can be deter¬ 
mined by repetitive computation of eigenvalues of a dynamic 
power grid model to give numerical sensitivities 0,0, ini, 
M- Also, there are exact computations of the sensitivity of 
the damping from a dynamic power grid model Q, lfl9l . 

ETIl that are based on the eigenvalue sensitivity formula 

dX _ wJpV 
op wv 

where w and v are left and right eigenvectors associated with 
the eigenvalue A and Jp is the derivative of the Jacobian with 
respect to the amount of generator redispatch p. The calcula¬ 
tion of Jp involves the Hessian and the sensitivity of the oper¬ 
ating point to p. However, requiring a large scale power system 
dynamic model poses difficulties. It is challenging to obtain 
validated models of generator dynamics over a wide area and 
particularly difficult to determine dynamic load models that 
would be applicable online when poor modal damping arises. 

We think that a good way to solve the difficulties with 
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online large scale power system dynamic models is to combine 
models with synchrophasor measurements to get actionable 
information about mode damping. In particular, the dynamic 
information about the power system can be estimated from 
synchrophasors. However, formula ([T]) is not suitable for this 
purpose since it is not feasible to estimate from measurements 
the left eigenvector w in 0 (or derivatives of eigenvectors 
in other versions of 0 )- This was a primary motivation for 
developing our new formula in im. In particular, the formula 
shows that the first order effect of a generator redispatch 
largely depends on the mode shape (the right eigenvector) 
and power flow quantities that can be measured online. The 
assumed equivalent generator dynamics only appears as a 
factor common to all redispatches. Given a lightly damped 
interarea mode of a system, the method can rank all the 
possible generators pairs. The rank is based on the size of 
the change of the damping ratio of the interarea mode. 

The main objective of ranking the generator pairs is to 
provide advice to the operator of several effective generator 
redispatches to damp the oscillations from which a corrective 
action can be selected. There are many economic goals and 
operational constraints governing the final selection of an 
appropriate dispatch by the operator, and the integration of 
this decision with optimal power flow and the power markets 
is left to future work. 

In this paper, the role in the method of measurements and 
the formula are first illustrated in a simple 3-bus system. Then 
the general formula derived in im is presented. How static 
and dynamic power quantities are obtained from measurements 
is discussed. The complex denominator of the formula contains 
the effect of the assumed equivalent generator dynamics, and 
a technique to obtain from measurements the phase of this 
complex denominator is presented. The paper then explains 
how to calculate the best generator pairs to damp a given 
mode, and illustrates and verifies the calculation with interarea 
modes of the New England 10-generator system. 

II. Formula for eigenvalue sensitivity with 

RESPECT TO GENERATOR REDISPATCH 
A. Special case of a 3-bus system 

We first present the eigenvalue sensitivity formula for a 
special case in which the formula simplifies. Consider the 
interarea mode of the simple 3-bus, 3-generator system shown 
in Fig. The generator dynamics is described by the swing 
equation and the transmission lines are lossless. In this special 
case, the bus voltage magnitudes are assumed constant. At the 
base case, the system is at an stable operating point and Fig.[T] 
shows the line power flows and the interarea oscillating mode 
pattern. The mode pattern shows that generator G1 is swinging 
against G3, and that G2 is not participating in the oscillation. 

The bus voltage phasor angles are i5i, (52, Let 6 i = S 1 —S 2 
and 02 = S 2 — S 3 be the voltage angle differences across line 1 
and line 2, and let pi and p 2 be the real power flows through 
line 1 and line 2. The interarea mode has complex eigenvalue A 
and complex right eigenvector x. Generator redispatch causes 
changes in the angles across the lines d 6 i and d 02 , and these in 


^-© 01 * • ©02 '■ ©G^ 

Fig. 1. The gray lines joining the buses show the magnitude of the 
power flow with the grayscale and the direction of the power flow 
with the arrows. The red arrows at each bus show the oscillation mode 
shape; that is, the magnitude and direction of the complex entries of 
the right eigenvector x associated with the interarea mode. 

turn cause changes dX in the eigenvalue. According to ifTTl . 
the formula for the first-order change in A with respect to 
generator redispatch in a 3-bus system with constant voltage 
magnitudes reduces to 

^ (x'ejVd6>i -f {x'g^fp^d 62 

a 

All quantities are in per unit unless otherwise indicated. The 
complex denominator a — 0.3080/87.12° seconds depends 
on the inertia and damping coefficients of the generators, the 
eigenvalue A and its right eigenvector x. The right eigenvector 
written with components corresponding to the bus angles is 
X ) The right eigenvector may also be 

expressed in terms of the angles across the line as x' = 
{x’g^.x’g^y. That is, x’g^ = xs^ — xs^ is the change in the 
eigenvector across line 1, and x'g_^ = xs^ — xs^ is the change 
of the eigenvector across line 2. It is these changes x'g^ and 
x'g_^ in the right eigenvector across the lines that appear in (j^. 

Formula ([^ is in terms of static load flow quantities 0 and p 
that are available from state estimation, and dynamic quantities 
A and x that could be available from synchrophasor measure¬ 
ments. Table |I] shows the values at the base case. Formula (|^ 
indicates which lines have suitable power flow and eigenvector 
components to affect oscillation damping. In particular, it is 
effective for the redispatch to change the angles across the 
lines that have both changes in the mode shape across the line 
and sufficient real power flow in the right direction. 


TABLE I 

Quantities from measurements 


Static quantities from 
state estimator 

Line No. 9 p [pu] 

Dynamic quantities 
from synchrophasors 

1 5.164° 0.2 

2 2.579° 0.1 

-0.2958 + j0.009838 
-0.2668 + j0.004026 


Let us define the complex coefficients of d0\ and d02 in 0 
as Cei Cg^ respectively. Then 

d\ = Cg^d0i Cg^d02, where (3) 

Ce^ = -0.000928 - j0.0569, Cg^ = 0.000462 - j0.0231. 

Coefficient Cg^ has the largest real and imaginary components. 
It follows that dX is more sensitive to redispatches done 
through line 1 than redispatches done through line 2. The 

'since x is the right eigenvector of a quadratic formulation of the 
eigenvalue problem EJ. it contains the angle components, but not the 
frequency components, of the conventional right eigenvector. That is, 
the conventional right eigenvector is = 

= (x,\xY. 
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damping ratio at the base case is given by 


Re{A} tr a 

|A| \a +ja;| 


(4) 


Table [I^ shows the results for the three possible generators 
pairs in the system for a small redispatch of 0.01 pu. As 
expected from ([^, the damping ratio has the largest changes 
when redispatch is done through line 1. This is the case of 
generator pairs G1+,G2- and G1+,G3-; their changes in d( 
are close. The damping ratio has almost no change when 
redispatch is done only through line 2 with the pair G3+,G2-. 
Gi+,Gj- indicates the redispatch that increases power at gen¬ 
erator Gi and decreases power at generator Gj. 


TABLE II 

Generator pairs ranked by change in ((%); redispatch = O.OI 

Generator Pair dQ{%) 

1 G1+,G2- 0.000393 

2 Gl-i-,G3- 0.000386 

3 G3-i-,G2- 0.000007 


B. General case of formula for the change in the eigenvalue 

Consider a connected power grid that has m generators, 
n + m buses and £ lines. Buses n + 1,..., n + m are the 
internal buses of the generators. Assume AC power flow 
and lossless transmission lines. Every generator is modeled 
with an equivalent second order equation (swing equation) 
and constant internal voltage magnitude. Loads are constant 
power, but frequency dependence of the real power and voltage 
dependence of the reactive power can be accommodated 
HD. Then the power system dynamics is described by a 
set of differential-algebraic equations with variables {5, V). 
The dynamic variables are ..., 6 n+m, and the algebraic 
variables are 6 i,... ,Sn,Vi,... ,Vn- In this paper we use the 
quadratic formulation of the eigenvalue problem HD, CD, 
||2^ with right eigenvector x and eigenvalue A: 

(MA^ + DX + L)x = 0, (5) 


where L is part of the system Jacobian described 
in HD and M and D are diagonal matrices con¬ 
taining the generator equivalent inertias and dampings; 
M = diag{2/ii/wo,2/i2,/a;o,...,2/im/wo,0, ...,0} and 
D = diagjdi, ^ 2 ,..., dm, 0, ■ • •, 0}. The components 
of the right eigenvector or mode shape are x = 

In other papers, the eigenstructure of differential-algebraic 
models with extended Jacobian J is often analyzed using 
the generalized eigenvalue problem Jv = XEv ll24ll . The 
difference between x and v is that x does not include the 
components of v related to generator angular speeds. 

According to HD, the new formula for the sensitivity of a 
nonresonant (algebraic multiplicity one) eigenvalue A is 



a 


Y, CeJOk + Y 


( 6 ) 


The denominator of (|^, which is the same for all redispatches 
of a given mode, is 

a = 2Ax^Mx -I- x^Dx, (7) 

dOk is the change in angle across the line k and dVi is the 
change in load i voltage magnitude due to the redispatch. 

C9k = + 2x'g^x'^^q^, ( 8 ) 


Gy, = E 

/c=l 


^7. = -L. 








Cp,P.)V-^ - CQ^QrV-\ (9) 



Pk = bkViVj sin 0k, 


qk = -bkViVj cos 0k, 


where 


bk 


5i — Sj if bus i is sending end of line k 
6j — Si if bus i is receiving end of line k 

Xj — x^ if bus i is sending end of line k 

Xj — x^ if bus i is receiving end of line k 

Xy 

—- H- - if line k joins load bus i to load bus j 

Xy. 

-f- if line k joins load bus i to generator bus j 


Pk is the real power flow in line k, qk is part of the reactive 
power flow in line k, and \ Aik \ is the absolute value of the ik- 
component of the bus-line incidence matrix. Qi is the reactive 
power demanded by load bus i. 

Formula ([^ expresses the first-order change in the eigen¬ 
value dX in terms of the changes d0 in angles across the lines 
and changes dV in load voltage magnitudes caused by the 
generator redispatch. dX depends linearly on the active power 
flow, part of the reactive power flow through every line of 
the network, and the reactive power demands at the loads, all 
evaluated at the base case. 


C. Generator modeling 

The overall dynamics of each generator is described by an 
equivalent swing equation model with constant internal voltage 
magnitude. For generator bus i, an “equivalent swing equa¬ 
tion” is the standard swing equation model with inertia and 
damping coefficients hi and di that produce the second order 
model that best approximates the entire dynamics of the gener¬ 
ator i and its controls (this would be described as a second or¬ 
der dominant pole approximation in automatic controls). In our 
approach, there is no need to model the parameters hi and di 
for each generator; it is only necessary to assume the existence 
of a second order model that can approximate well enough the 
contribution of the generator to the electromechanical mode. 
Indeed, since (j^ only includes the combined generator dynam¬ 
ics as a common factor that is the same for all redispatches, 
we do not need to know the individual parameters of each 
equivalent generator model. Other authors using synchropha- 
sor measurements to identify aggregated generator dynamics 
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for studying oscillations also assume swing equation generator 
models but for their purposes identify the individual inertia and 
damping coefficients ll25l . Ebll . However, we have to consider 
the generator modeling independently of previous studies since 
our generator redispatch application has novel and different 


modeling requirements as discussed in section VII 


III. Measurements 

Formula (|^ depends on power system quantities that could 
be observed from measurements from the state estimator and 
from synchrophasors. 


A. Quantities obtained from state estimator: p, q, Q, d9, dV 

The state estimator can determine the active power p and 
part of the reactive power flow q through every line of the 
network, and the load reactive power injections Q. Load flow 
equations can be used to relate the generator redispatches 
to changes in angles across the lines dO and load voltage 
magnitudes dV. 

B. Quantities obtained from synchrophasors: \,x 

Formula 0 depends on the dynamic quantities A and 
X that satisfy the quadratic formulation of the eigenvalue 
problem Q. For an electromechanical oscillation present in 
a system, synchrophasors can make online measurements 0, 
Gol, Ell, ms) of the damping and frequency of the eigenvalue 
A associated with the oscillation, x is easy to obtain from a 
conventional right eigenvector. The right eigenvector of A is in 
principle, and to some considerable extent in practice, avail¬ 
able from ambient or transient synchrophasor measurements 
El, ca, Qa. It is conceivable that historical observations 
or offline computations or general principles about the mode 
shape could be used to augment or interpolate the real-time 
observations, or that the real-time observations could be used 
to verify a predicted mode shape. Thus some combination of 
measurements and calculation could yield the mode shape x. 


C. Method for estimating phase of a from measurements 

The complex denominator a = 2Xx'^Mx + x^Dx is 
the only part of formula 0 that depends on the generator 
equivalent dynamic parameters in the inertia and damping 
matrices M and D. In estimating the change dX in the 
eigenvalue from (|^, the angle Za controls the direction of 
dX and the size of a controls the size of dX. For each mode, 
a is the same for all generator redispatches. Therefore, in 
order to rank the generator redispatches for a given mode, 
it is sufficient to estimate Za from measurements. 

Formula 0 can be summarized as 

Numerator 

dX = -, so Za = ZNumerator — ZaA. (10) 

a 

To estimate Za we propose to take advantage of the small 
random load variations around the operating equilibrium. 
For such small random load variations, samples of dX could 
be obtained from a series of synchrophasor estimates of A. 
Samples of dO and dV can be obtained from the load flow 


equations with simulated random load variations, so that 
samples of the numerator of 0 can be computed. It turns out 
that both the dX samples and numerator samples have preferred 
or principal directions in the complex plane. Analyzing the dX 
and numerator samples with Principal Component Analysis 
gives a principal axis direction for dX and a principal axis 
direction for the numerator, and according to ( [T0| ), the angle 
between these principal axis directions can be used to find 
Za. In section VI-B[ we illustrate and apply this calculation 
of Za to interarea modes of the New England system. 


D. Measurement processing requirements 

To be able to supply the dynamic quantities in formula 
0, online synchrophasor monitoring of the oscillatory mode 
eigenvalue A and mode shape x is needed. The modal eigen¬ 
value is used directly in 0 and, according to the suggested 
approach in section |III-C| for estimating the phase of a. This 
subsection comments briefly about the likely measurement 
processing limitations. We would expect to be able to use 
standard sampling rates and data concentration similar to that 
already in use for mode monitoring. 

Methods to estimate oscillatory modes and mode shapes 
from synchrophasor data are deployed and improving, and 
some recent methods ca, m, Eoi, ED, have used up to 
5 minutes of ambient or transient data for multiple parallel 
algorithms to converge to consistent results. There is also con¬ 
sistency between results from methods for transient response 
after a disturbance and ambient methods. However, it does not 
matter for our algorithm whether the dynamics is estimated 
from transient or ambient responses. 

The mode shape estimation currently samples the mode 
shape at multiple spatial locations. We require the mode 
shape at other locations to be interpolated, perhaps guided 
by previously observed mode shapes for specific modes. Our 
approach also requires static quantities to be estimated with 
standard state estimation. Given state estimator convergence, 
the estimated state should be readily available within the 
several minute time scale already required for mode estima¬ 
tion. We do not anticipate significant computational delay in 
evaluating formula (j^ and ranking the generators. 


IV. Relating redispatch to the changes dO, dV, dX 
Formula (j^ expresses the eigenvalue change dX in terms of 
dO and dV. It remains to express dd and dV in terms of the 
redispatch dP. Define the coefficient vectors of dO and dV in 
(|^ as 

Cg = —Cg and Cy = —Cy. (11) 

a a 

Then 

dX = Cg-dO + Cy -dV = (Ce, Cy) . (12) 

From im we know that the linear relationship between 
the changes in angles d9 and changes in voltages dV with 
redispatch dP is given by 


( fdP 

\dV) yO„x(n-|-m) InxnJ \ ^ , 


( 13 ) 
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where A is the network incidence matrix, and f indicates 
pseudoinverse. As an alternative to the linearized computation 


of {dO,dV) from the redispatch dP in (13i, one can simply 


recompute the AC load flow with the assumed change in 
redispatch dP. We can relate the change in an eigenvalue of 


a mode to a generator redispatch dP by substituting (131 into 

([T2J: 


dX = Cp-dP = J2 Cp^dP,. 


(14) 


The complex coefficient Cp. gives the contribution of gener¬ 
ator i to dX by increasing the real power of generator i by 
dPi. The redispatch is assumed to satisfy the active power 
balance constraint dPi = 0. It is convenient to define 
Cp = |a|Cp. Then multiplying both sides of (14i by |a| 
gives 


\a\dX = Cp-dP^y"Cf 


dP. 


(15) 


It can be seen by considering (lli-(14i that Cp can be 
calculated from Cg, Cy (Note, for example, that 

\a\Cg = e~^°‘Cg.) Since |a|dA is dX multiplied by a real 
constant, the complex number |a|dA has the same angle as 
dX and has magnitude proportional to dX. It follows that we 
can use jajdA to rank the generator redispatches. 


V. Ranking generators pairs by increase in 

DAMPING RATIO DUE TO REDISPATCH 


For a given oscillatory mode with insufficient damping ratio, 
the first step towards ranking the best generators to redispatch 
is to make the following measurements and calculations; 

1) The eigenvalue A associated with the interarea oscil¬ 
lation and its mode shape x are estimated from syn- 
chrophasor measurements. 

2) p, q, Q are obtained from the state estimator. 

3) Za is computed from measurements. 

4) Coefficients Cg^ and Cy. in (|^ are computed from ^ 
and 

5) Cp is computed from Cg, Cy, Za. 


The most straightforward way of implementing generator 
redispatch is with a pair GiH-,Gj- of generators; that is, in¬ 
creasing the real power of generator i by some amount dPi 
and increasing the real power of generator j by the negative 
amount dPj = —dPi. We can use (141 to find the generator 
pair Gi-i-,Gj- that gives the most favorable eigenvalue change 
dX that best increases the damping ratio. For a redispatch 
of amount dPi with generator pair Gi-i-,Gj-, the resulting 
eigenvalue change satisfies 


\a\dKj=[Cp,-Cp^)dP,. (16) 


Also the opposite redispatch Gj-,Gi-i- with the same generator 
pair will yield \a\dXji = —\a\dXij. For every pair of gener¬ 
ators in the system we can compute \a\dXij and \a\dXji and 
then we can compute the corresponding damping ratios Qj 
and Qi for every generator pair with (|^. The damping ratios 
for all pairs {C 12 , C 21 , • ■ ■, Cm.m-i} are then ranked 


from the largest one to the smallest one. The highly ranked 
pairs indicate the generator pairs and redispatch directions that 
will increase the damping ratio of the mode the most. 


VI. Damping modes of the New England system 

This section presents the results for damping interarea 
modes of the New England 10-generator system 0^ with 
generator redispatch. The equivalent parameters of the gen¬ 
erators are given in Appendix A and the rest of the data is 
provided in lf3^ . All the numerical computation was done with 


the software Mathematica. Table III shows the eigenvalues 


of the 4 interarea modes at the base case. Since the system 
has 10 generators, there are 45 generators pairs. The method 
presented in section was implemented for the interarea 
modes (this computation used a calculated value of a and 
a recalculated load flow to evaluate d9 and dV from dP). 


TABLE III 

Interarea mode eigenvalues at the base case 


Mode Eigenvalue A [1/s] f [Hz] Damping Ratio ((%) 

1 -0.040336 + j3.4135 0.54327 1.18157 

2 -0.018839 -I- j4.7631 0.75807 0.39551 

3 -0.024903 -I- j5.4994 0.87526 0.45283 

4 -0.055799 -i- j6.0159 0.95746_0.92748 


A. Verifying the formula 

Eormula was verified for the four interarea modes, but 
here we present results only for mode 1. The eigenvalue Ai 
was computed for a given redispatch using (j^ and using the 
exact computation of the Jacobian at the new operating point. 
Table IV shows the exact and approximate Ai for different 
amounts of redispatch of generator pair G5-i-,G9-0 

Table |IV] confirms that (|^ reproduces the first order varia¬ 
tion of Ai with respect to the redispatch. 


TABLE IV 

Eigenvalue Ai for generator redispatch G5-I-, G9- 


Redispatch 

Exact eigenvalue 

Approximate eigenvalue 

0.0 

-0.0403355 - 1 -j3.4135 

-0.0403355 -H j3.4135 

0.0005 

-0.0403225 - 1 -j3.4131 

-0.0403225 -H j3.4131 

0.001 

-0.0403095 - 1 -j3.4127 

-0.0403095 + j3.4127 

0.01 

-0.0400715 - 1 -j3.4059 

-0.0400736 + j3.4059 

0.02 

-0.0397985 - 1 -j3.3981 

-0.0398086 + j3.3983 

0.03 

-0.0395161 -i-j3.3899 

-0.0395403 + j3.3905 


B. Estimating the phase of a for random load variations 

Table |V] shows the exact Zct computed with (0 from the 
mode shape and the equivalent generator parameters, and the 
Za estimated from random load variations with the method 
of section ITlI-CI 

^ The base case generations of the generators G1 through GIO are 
2.5,4.8,6.5,6.32,5.08,6.5,5.6,5.4,8.3,10.0 per unit respectively. As an 
example of specifying the redispatch, a 0.03 redispatch of generator pair 
G5 -h,G 9- changes the generation of G5 from 5.08 to 5.11 and the generation 
of G9 from 8.3 to 8.27. 
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For each interarea mode, the set of random loads used in the 
computation were generated with the software Mathematica. 
The active power random load vector P'' was sampled from 
a normal distribution of zero mear0 The components of the 
reactive power random load vector Q'^ were computed as 

q : = (17) 

Pi and Qi are the active and reactive power demanded by 
load i at the base case. Then the load flow solution for the 
vector of loads P + P^ and Q + Q'' were computed and 
d\ and the numerator of were computed. 50 random 
load scenarios were generated. Fig. shows samples for Ai 
after being trimmed by 30% to remove outliers and analyzed 
with principal component analysis. This 2-dimensional Ai 
data was trimmed with multidimensional trimming based on 
projection depth 1^ . Projection depth induces order for high 
dimensional data, which makes trimming straightforward. The 
results in Table |V] show that the method gives a very good 
estimation of Za. 


TABLE V 

Estimated and exact phases of a for the interarea modes 


Mode 

Exact ZcK 

Estimated Za 

Exact Zq! - Estimated Za 

1 

88.658° 

88.718° 

-0.0600° 

2 

90.257° 

90.281° 

-0.0241° 

3 

89.700° 

89.685° 

0.0156° 

4 

89.026° 

89.126° 

-0.1004° 



Numerator 



Re{dA,) 


Eig. 2. 50 samples of the numerator of and dX after trimming 

by 30%. Principal axes are computed and shown as lines. 


C. Ranking generator pairs to damp mode 1 

Fig. shows the power flow p and oscillating mode pattern 
of interarea mode 1 at the base case. The mode pattern shows 
that generators G2 through G9 are oscillating against GIO. 
The component of GIO is not very large compared with 
the components of the other generators, but GIO is a large 
generator that represents an equivalent of the New York State 
grid. From Fig. we can see that generator G5 participates 
most in the oscillation, G8 has a small participation in the 
oscillation, and G1 does not participate in the oscillation. The 
method of section was applied to rank the generator pairs 
that best increase the damping ratio for a small redispatch of 


0.01 pu, and the top 10 generator pairs are shown in Table VI 


The top 9 generator pairs all include G5-; that is, decreasing 
generation at G5 with increasing generation elsewhere. The 
change in damping ratio for these pairs are of the same order of 
magnitude, so it is clear from Table |VI| that the largest changes 
in damping ratio are due to the generator pairs involving G5-. 


TABLE VI 

Generator pairs ranked by change in Ci(%); redispatch = 0.01 



Eig. 3. The gray lines joining the buses show the magnitude of 
the active power flow with the grayscale and the direction of the 
active power flow with the arrows. The red arrows at each bus show 
the pattern of the oscillation for mode 1; that is, the magnitude and 
direction of the entries of the right eigenvector xs- 


Generator Pair 

dCl(%) 

Generator Pair 

dCi(%) 

1 G6+,G5- 

0.00554 

6 G1+,G5- 

0.00505 

2 G3+,G5- 

0.00549 

7 G9+,G5- 

0.00502 

3 G7+,G5- 

0.00545 

8 G4+,G5- 

0.00489 

4 G2+,G5- 

0.00537 

9 G10+,G5- 

0.00481 

5 G8+,G5- 

0.00510 

10 G6+,G10- 

0.000725 


^Load 12 has exceptionally high reactive power and was not vaiied. 


D. Ranking generator pairs to damp mode 2 

Fig. 1^ shows the power flow p and oscillating mode pattern 
of interarea mode A 2 at the base case. The mode pattern shows 
that generators G1-G3 and G6-G9 are oscillating against G5. 
Generators G5 and G9 participate most in the oscillation. 

The method of section |V] was applied to rank the generator 
pairs that best increase the damping ratio for a small redispatch 
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Fig. 4. The gray lines joining the buses show the magnitude of the 
active power flow with the grayscale and the direction of the active 
power flow with the arrows. The red arrows at each bus show the 
pattern of the oscillation for mode 2. 


of 0.01 pu, and the top 10 generator pairs are shown in 


Table VII The top 9 generator pairs all include G5+; that 
is, increasing generation at G5 with decreasing generation 
elsewhere. Although the change in damping ratio is of the 
same order of magnitude for the 10 pairs, the increase in 
damping ratio for pair 10 G4+,G9-, is roughly half of the 
increase of damping ratio for pair 9 G5+,G4-. 

Now we analyze the most significant components of formula 
^ to explain why G5 is playing a key role in damping mode 2. 
Table |VII| shows that the pairs with the largest change in 


TABLE VII 

Generator pairs ranked by change in C2(%); redispatch = 0.01 


Re{Cv} 


Gen. Pair 

dC,2(%) 

dX 

Re{Ce} ■ de 

■dV 

1 G5+,G9- 

9.60E-3 

-4.23E-4 - jO.008 

-3.04E-4 

-1.19E-4 

2 G5+,G8- 

7.30E-3 

-3.06E-4 - jO.OlO 

-2.23E-4 

-8.32E-5 

3 G5+,G1- 

7.28E-3 

-3.02E-4-j0.011 

-2.21E-4 

-8.12E-5 

4 G5+,G10- 

7.26E-3 

-3.00E-4-j0.011 

-2.19E-4 

-8.10E-5 

5 G5+,G2- 

6.85E-3 

-2.85E-4 - jO.OlO 

-2.10E-4 

-7.55E-5 

6 G5+,G7- 

6.79E-3 

-2.84E-4 - jO.OlO 

-2.11E-4 

-7.30E-5 

7 G5+,G6- 

6.77E-3 

-2.84E-4 - jO.OlO 

-2.11E-4 

-7.28E-5 

8 G5+,G3- 

6.74E-3 

-2.81E-4 - jO.OlO 

-2.08E-4 

-7.38E-5 

9 G5+,G4- 

6.05E-3 

-2.52E-4 - j0.009 

-1.90E-4 

-6.21E-5 

10 G4+,G9- 

3.53E-3 

-1.71E-4 - jO.OOl 

-1.14E-4 

-5.66E-5 


damping ratio dC are the ones with the largest increase in the 
damping Re{dA}, so we can focus on the real part of (12i. 
Moreover, Table VII also shows that the changes RejCg} • dO 
are larger than the changes Re{Cy} • dV, so we focus on 
analyzing the terms of related to d9: 

Re{dA} = RejCe} • d9 + Re{Cy} • dV « RejCe} • d9 (18) 


Fig. 0 shows different quantities related to the 56 lines of 
the New England system in gray scale. Each generator is 
represented in the network by an internal bus and a terminal 
bus. There are two lines in series at the edge of the network 
associated with each generator. The line joining the internal 
bus to the terminal bus represents the generator transient 


o o 




(c) \de\ for G5+,G9-. 


(d) |Re{Ce}-d6»| for A 2 & G5+,G9-. 


Fig. 5. Gray scale in the lines shows the components of the specified vector. 


reactance, and the other line lumps together the transformer 
and lines joining the generator terminal bus to the network. 
(GIO differs since it represents New York state.) Eig. 5|a ) 
shows iRelCell, the absolute value of d9’s coefficient in (18 1 . 
The lines that represents the transient reactance of G5 and the 
transient reactance of G9 have large components of iRelCg}!. 
Eig.|^b) shows \dp\, the absolute value of the change in power 
flow in lines for the best ranked generator pair G5+,G9-. As 
expected, several lines have a significant change in power 
flow, but Eig. I^c) shows that only the line that represents 
the transient reactance of G5 has a large change in angle 
\d9\. This is due to the fact that transient reactance of G5 
is much larger than the reactance of any of the other 55 
lines of the system (see appendix A). So, for any possible 
generator pair that involves G5, the line that represents the 
transient reactance of G5 will always have the largest change 
in angle. This large change in angle, combined with a large 
coefficient iRelCe}!, produces the dominant term of ( [T^ , as 
shown by Eig.|^d). Thus G5 is the key generator to participate 
in redispatch for producing the largest changes in damping for 
mode 2. Although Eig. [^a) shows that there are other lines 
that have large d9’s real coefficient, and Eig. |^b) shows that 
other lines have an important change in power \dp\, Eig. |^c) 
shows that such lines do not have a large change in angle d9, 
and as a result their associated terms in ( fTS] ! are not large, as 
seen in Fig. [^d). This analysis shows a mechanism of how 
damping by redispatch works by changing the angles across 
lines that have large coefficients Cg. 












E. Larger redispatches 

Larger redispatches introduce some nonlinearity in the 
change in the damping ratio. This nonlinearity arises from 
three sources: the change in the load flow, the change in 
the eigenvalue dX, and the change in the damping ratio. 
For a given size of redispatch of a generator pair, we can 
compute the exact nonlinear change in the load flow due 
to the redispatch, use (|^ to linearly estimate the change in 
the eigenvalue based on measurements, and then compute the 
nonlinear change in the damping ratio. This is a large signal 
application of (|^, and the top 10 generator pairs for modes 1 
and 2 are shown in Table |VIII| for a redispatch of 0.4 pu. For 
comparison. Table |VIII also shows the exact calculation that 
uses the nonlinear computation of dX. As shown in Fig. the 
use of (|^ for a larger redispatch gives the same grouping of 
the top 9 effective generator pairs and a similar ranking and 
grouping of the top 10 generator pairs as the exact calculation. 
An exception is that for mode 1, the 5th ranked generator pair 
G9+,G5- for a larger redispatch using becomes the 9th 
ranked generator pair for the exact calculation. 

Fig. i also compares the changes in the damping ratio of 
the top 10 generator pairs for a small redispatch using (6) 
to the larger redispatch using (6). The grouping of the top 
9 effective generator pairs is preserved and the approximate 
ranking is preserved. 

Overall, some details of the rankings differ for similarly 
effective generator pairs, but since the ranking will be used 
to provide a set of effective generators pairs from which an 
operationally suitable pair can be selected for redispatch by 
operators, the performance of the ranking is satisfactory. 


TABLE VIII 

Generator pairs ranked by change in C(%); redispatch = 0.4 


MODE 1 
Gen. Pair 

dG 

Using 

(%) 

Exact 

MODE 2 
Gen. Pair 

dG 

Using 

(%) 

Exact 

1 G6+,G5- 

o.isP 

0.138 

1 G5+,G9- 

0.563 

0.600 

2 G7+,G5- 

0.177 

0.134 

2 G5+,G10- 

0.558 

0.615 

3 G3+,G5- 

0.177 

0.135 

3 G5+,G1- 

0.534 

0.595 

4 G2+,G5- 

0.171 

0.128 

4 G5+,G8- 

0.506 

0.572 

5 G9+,G5- 

0.162 

0.092 

5 G5+,G2- 

0.473 

0.545 

6 G4+,G5- 

0.160 

0.122 

6 G5+,G3- 

0.453 

0.527 

7 G8+,G5- 

0.159 

0.115 

7 G5+,G7- 

0.434 

0.506 

8 GI+,G5- 

0.155 

0.112 

8 G5+,G6- 

0.432 

0.503 

9 GI0+,G5- 

0.145 

0.104 

9 G5+,G4- 

0.365 

0.444 

10 G6+,GI0- 

0.030 

0.034 

10 G4+,G9- 

0.139 

0.120 


VII. Discussion of generator modeling 


Since our approach depends in new ways on both mea¬ 
surements and an equivalent second-order swing equation 
generator model, our generator modeling requirements are 
different than in other approaches to suppressing oscillations. 


Section II-C explains that the generator dynamics are approx¬ 
imated by a swing equation, but we do not need to determine 
the parameters of the swing equation for any individual gen¬ 
erator. This section further discusses the generator modeling. 

As a general observation, in closed loop control of oscil¬ 
lations with power system stabilizers, which forms most of 
the literature on suppressing oscillations, the generator and its 
controls need to be modeled in sufficient detail. Indeed the 


Using formula (6) 


MODE 2 

Using formula (6) 

A Exact 



Redispatch 0.01 pu 


-Y- 

Redispatch 0.4 pu 


Redispatch 0.01 pu 


Redispatch 0.4 pu 


Fig. 6. Changes in damping ratio for Mode 1 and Mode 2 for the top 
10 generator pairs from Tables [VT] |VII| [VIII| rescaled to the same range to 
allow comparison. For each mode, left hand dots use formula 0 for a small 
redispatch, middle dots use for a larger redispatch, and the right hand 
dots are the exact calculation for the larger redispatch. Clustering of similarly 
effective generator pairs is shown by close dots and changes in ranking appear 
as lines crossing. 


designed control gains directly affect entries of the Jacobian 
to damp the oscillatory mode. Our control is open loop and 
works by the entirely different principle of exploiting system 
nonlinearity by changing the operating point at which the 
Jacobian is evaluated. This changes the focus from the linear 
parts of the model to the nonlinearities. 

Formula (j^ computes the first order sensitivity of the oscil¬ 
latory mode eigenvalue to generator redispatch. Appendix B 
proves that the first order eigenvalue sensitivity to redispatch 
does not depend on linear parts of the power system model. In 
particular, if the generator magnetic saturation and hysteresis 
are neglected, the higher order parts of the generator modeling 
are linear and the eigenvalue sensitivity only depends on the 
nonlinearity in the swing equation and any stator algebraic 
equations and does not depend on the linear higher order part 
of the generator dynamic modeling. This does suggest that 
the linear higher order generator dynamics can be omitted in 
deriving the formula. 

In applying formula we do not use a model of the 
power system dynamics. Instead we rely on measurements 
of the system dynamics, particularly the eigenvalue and the 
right eigenvector of the mode and the phase of the complex 
scalar parameter a that combines together all of the generator 
dynamics. There are no model assumptions in these measured 
quantities. That is, if part of the power system affects the oscil¬ 
latory dynamics, the effect will appear via the measurements 
used by the formula. 

We did an initial test of the approximation involved in 
generator modeling with the interarea mode of the 3-generator 
model similar to Fig. 1 with a sixth-order generator model at 
each bus. Formula is applied to this detailed model by 
measuring the phase of a, and using the computed mode, 
mode shape, and load flow to estimate the change dX in 
the eigenvalue that would arise from small redispatches in 
generation. Then the detailed model is used to compute the 
exact change dX in the eigenvalue. The comparison of the 
approximate and exact dX is shown in Table IX Similarly 
to the intended application of the formula to ranking redis¬ 
patches in a real power system, the power system model 
with sixth-order generators does not have the parameters of 
the equivalent second order generator models available. That 
is, the magnitude of a is not known, and only the phase of 
a is estimated, and so the formula predicts dX to within a 
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constant real multiplier. Therefore in Table we compare 
the ratios of |c?A| for each redispatch to |dA| for the redispatch 
of generators G1 and G2, as well as comparing the phases of 
d\. The approximation of d\ in Table is close enough to 
be acceptable for ranking of generator redispatches. 

While the generator modeling issues should be investigated 
further in future work, and further analytic progress is not 
ruled out, both theoretical considerations of the irrelevance of 
linear parts of the generator model and an initial test indicate 
that combining measurements of the dynamic quantities with 
a formula assuming a second order swing equation can be 
adequate for ranking generator redispatches. 

TABLE IX 

Eigenvalue changes for redispatch of 0.01 pu in 3-generator 

SYSTEM WITH SIXTH-ORDER GENERATOR MODELING 


Generator 

pair 

Exact 

ZdA \dX\ ratio 

Approximate with formula 
ZdX |dA| ratio 

G1+,G2- 

-109.5° 

1.00 

-106.3° 

1.00 

G2+,G3- 

-98.6° 

2.54 

-98.5° 

2.12 

G1+,G3- 

-101.7° 

3.53 

-101.0° 

3.11 


VIII 

Conclusions 



There has been success in monitoring interarea modal 
damping with synchrophasor measurements El, Col; the next 
step is to leverage synchrophasor measurements to provide 
advice to the operators to maintain a suitable modal damping 
ratio when the damping ratio is insufficient. Difficulties in 
accomplishing this in the past include the lack of wide-area 
online dynamic models and standard formulas that depend 
on quantities that cannot be measured. In this paper, we 
circumvent these difficulties by calculating the best generator 
pairs to redispatch to maintain modal damping by combining 
synchrophasor and state estimator measurements with a new 
analytic formula for the sensitivity of the mode eigenvalue with 
respect to generator redispatch. The assumed equivalent gener¬ 
ator dynamics only appears as a complex factor 1 /a common 
to all redispatches and we propose a method of estimating the 
phase of a from ambient measurements. The new formula is 
somewhat complicated, and we explain and illustrate how it 
works in 3 and 10 generator examples. Future work may well 
discover further insights and applications using the formula. 
In summary, we make substantial progress towards practical 
application of a new formula to damp interarea oscillations 
based on measurable quantities. 
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Appendix A: New England generator data 


Gen. 

No. 

Terminal 
Bus No. 

Internal 
Bus No. 

h [s] 

d [s] 

< 

V Internal 
Bus 

1 

30 

40 

42.0 

0.0267 

0.031 

1.0501 

2 

31 

41 

30.3 

0.0161 

0.0697 

1.0388 

3 

32 

42 

35.8 

0.0209 

0.0531 

1.0439 

4 

33 

43 

28.6 

0.0243 

0.0436 

1.0348 

5 

34 

44 

26.0 

0.0014 

0.132 

1.2098 

6 

35 

45 

34.8 

0.0277 

0.05 

1.0941 

7 

36 

46 

26.4 

0.0140 

0.049 

1.0944 

8 

37 

47 

24.3 

0.0116 

0.057 

1.0705 

9 

38 

48 

34.5 

0.0002 

0.057 

1.1252 

10 

39 

49 

500.0 

0.3979 

0.006 

1.0317 


Vi is the internal constant voltage magnitude of generator i at the base case; 
i.e., at zero redispatch, and is the transient generator reactance. 


Appendix B: Irrelevance of linear modeling 


This appendix proves from Q that the first order eigenvalue 
sensitivity to redispatch does not depend on linear parts of the 
power system model. This result was first mentioned in ||7] 
section 4.6]. It is convenient to use the extended differential- 
algebraic form of the power system equations ll24ll with state 
vector z, Jacobian J, and extended right and left eigenvectors 
V and w. In this notation, (1) becomes 


dX wJpV 
dp wv 


(19) 


In this case, the parameter p parameterizes the generator 
redispatch, and it appears linearly in the system equations. 
Therefore p does not appear explicitly in the Jacobian J, and 


7 = — = V 
p dp ^ dz>^ p 

k 


( 20 ) 


where Zp = is the sensitivity of the operating point z to 
the redispatch. Then substituting in ( [T9| ) and writing it out in 
coordinates gives 


dp 


k-i 




( 21 ) 


It is clear from (21 1 that the linear parts of the model vanish 
in and that if a dynamic state is associated with a 

linear differential equation, the corresponding entry of w gets 
multiplied by zero in the numerator of (|2l]i. 
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